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ABSTRACT 


In  order  to  validate  existing  techniques  for  computer 
simulation  of  the  welding  process,  an  experiment  is  designed 
and  carried  out  in  which  seam  and  girth  welds  are  produced 
on  a  high  strength  steel  cylinder.  Fabrication  of  the  HY-130 
cylinder,  selection  of  welding  parameters  for  the  seam  and 
girth  welds,  and  the  selection  and  arrangement  of  instrumenta¬ 
tion  are  described.  The  theoretical  background  of  finite 
element  computer  simulation  of  heat  flow  in  weldments  as  well 
as  techniques  for  axisymmetric  finite  element  analysis  of 
thermal  stress  and  metal  movement  during  welding,  are  dis¬ 
cussed.  Finite  element  representation  of  weldment  geometry 
for  heat  flow  and  stress  analysis  is  presented.  Metallurgical 
properties  of  the  HY-130  cylinder  seam  weld  are  presented 
in  the  form  of  photo-micrographs  and  electron  scan  micro¬ 
graphs  . 

Temperature,  strain,  and  radial  distortion  are  measured 
at  various  locations  on  the  HY-130  cylinder  and  the  measured 
values  are  compared  in  tabular  and  graphical  form  with  those 
values  predicted  by  use  of  the  heat  flow  and  stress  analysis 
programs.  Excellent  agreement  is  achieved  between  measured 
and  analytically  predicted  temperature  distributions,  but 
only  qualitative  agreement  is  shown  for  strains  and  radial 
distortion. 
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CHAPTER  I 

I. A  INTRODUCTION: 

Welding  has  been  used  as  the  primary  means  of  fabricating 
ocean  structures  for  the  past  forty  years.  Design  tools  for 
the  fabrication  of  welded  structures  have  evolved  primarily 
as  the  result  of  experimentation.  These  tools  are  usually 
expressed  as  empirical  relationships  among  welding  parameters, 
some  of  which  include  structure  and  weld  joint  geometry, 
physical  properties  of  the  structural  materials  and  welding 
consummables ,  as  well  as  methods  and  rates  of  energy  trans¬ 
fer  during  the  joining  process.  In  most  cases  the  use  of 
these  empirical  relationships  has  been  adequate  to  ensure  the 
achievement  of  good  quality  weld  joints  in  a  production  en¬ 
vironment. 

As  design  tools  become  more  sophisticated,  and  the  in¬ 
terrelationships  among  physical  properties  of  welding  and 
structural  materials  are  better  understood,  design  practices 
relating  to  allowable  stresses  in  welded  structures  converge 
toward  material  limits.  Development  and  refinement  of  tech¬ 
niques  for  detecting  and  classifying  flaws  in  weldments  is 
another  design  consideration  which  tends  to  increase  confi¬ 
dence  in  welded  structures  subjected  to  more  demanding  re¬ 
quirements,  allowing  design  engineers  to  increase  allowable 
stresses  or  to  reduce  factors  of  safety  or  uncertainty. 


An  area  of  welding  technology  which  requires  further 
definition  and  refinement  is  the  occurrence  of  residual 
stresses  and  structural  distortion  during  and  as  the  result 
of  the  welding  process.  Extreme  non-uniform  thermal  cycles 
of  weld  metal  and  of  structural  metal  in  areas  near  the 
weldment  produce  unevenly  distributed  thermal  stresses 
which  cause  localized  plastic  deformations/  with  resultant 
residual  stresses  and  structural  distortion.  These  residual 
stresses  can  significantly  contribute  to  structural  failure, 
particularly  when  the  failure  mechanism  is  dependent  upon 
low  applied  stress,  such  as  brittle  fracture  or  failure  re¬ 
sulting  from  stress  corrosion  cracking.  Distortions  caused 
by  welding  can  also  cause  structural  failure  by  reducing  a 
structure's  ability  to  withstand  buckling  loads. 

The  ability  to  accurately  predict  the  magnit.de  and 
type  of  residual  stresses  and  distortion  resulting  from  the 
welding  of  ocean  structures  would  greatly  enhance  a  design 
engineer's  ability  to  estimate  total  stresses  which  materials 
must  withstand  during  loading.  More  accurate  estimation  of 
total  stresses  would  allow  for  more  effective  use  of  materials 
and  the  increased  predictability  of  structural  failure  would 
ultimately  result  in  the  production  of  more  reliable  ocean 


structures . 
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I.B  BACKGROUND; 

Attempts  to  analytically  predict  or  to  simulate  thermal 
stresses  and  metal  movement  which  occur  during  the  welding 
process  have  been  the  subject  of  considerable  research  effort 
during  the  past  ten  years.  The  complexity  of  the  interrela¬ 
tionships  among  welding  parameters  which  contribute  to  the 
formation  of  residual  stresses  in  weldments  does  not  allow 
the  development  of  empirical  relationships  which  are  useful 
for  more  than  one  type  of  material,  joint  geometry,  or  welding 
process.  Welding  variables  were  too  numerous  and  their  in¬ 
terrelationships  too  complex  for  the  development  of  a  general 
technique  for  the  analysis  of  transient  thermal  stresses, 
residual  stresses,  and  distortion  in  weldments.  Design 
engineers  usually  relied  upon  trial  and  error,  or  the  use  of 
design  lanes  derived  from  similar  applications  of  proposed 
welding  processes. 

Computer  simulation  of  the  welding  process  has  offered 
reasonable  opportunities  to  develop  techniques  which  apply 
to  wide  ranges  of  materials,  weld  geometries,  and  processes. 

A  description  of  the  chronological  and  systematic  development 
of  techniques  used  for  analysis  of  thermal  stresses  and  metal 
movement  during  the  welding  process  is  contained  in  a  paper 
presented  by  Masubuchi^  to  an  international  conference  on 
Residual  Stresses  in  Welded  Construction  and  their  Effects. 


This  conference  was  held  in  London  15-17  November  1977, 
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hosted  by  the  Welding  Institute,  Abington  Hall,  Abington, 

Cambridge,  England.  Systematic  development  of  computer 

techniques  for  welding  simulation  began  with  a  one-dimen- 

2 

sional  analysis  by  Tall  in  1961.  This  research  involved 
the  simulation  of  thermal  stresses  occurring  during  the 
bead  welding  of  a  strip.  Efforts  at  Massachusetts  Institute 
of  Technology  by  Masubuchi^  in  1968  expanded  Tail's  analysis 
into  a  one-dimensional  FORTRAN  program.  This  research  was 
expanded  into  the  development  of  a  two-dimensional  FORTRAN 

4 

program  in  1970  by  Iwaki  .  Iwaki  attempted  to  apply  tv j- 

dimensional  finite  analysis  to  bead  on  plate  welding.  The 

most  recent  computer  program  to  be  used  at  MIT  was  developed 

by  Muraki5.  Muraki's  program  uses  finite  element  analysis 

• 

techniques  to  predict  thermal  stresses  and  metal  movement 
based  upon  temperature  distributions  within  the  weldment 
which  are  calculated  by  heat  flow  programs  also  developed 
at  MIT. 

I .C  SCOPE  AND  OBJECTIVES  OF  THIS  STUDY: 

The  purpose  of  this  research  is  to  design  and  carry 
out  an  experimental  procedure  to  generate  data  which  can 
be  used  to  validate  the  existing  heat  flow  program  and  the 
program  which  uses  finite  element  analysis  to  predict  the 
occurrence  of  thermal  stresses  and  distortion  during  welding. 
The  production  of  seam  and  girth  welds  in  a  high  strength 
steel  cylinder  under  very  carefully  controlled  conditions 


while  measuring  temperature,  strain,  and  distortion  as  a 
function  of  time  and  location  with  respect  to  arc  travel 
would  serve  as  a  valuable  aid  in  evaluating  the  accuracy 
and  potential  usefulness  of  the  heat  flow  program  and  of 
Muraki's  axisymmetric  finite  element  program  for  analysis 
of  thermal  stress  and  metal  movement  during  welding.  This 
research  is  part  of  a  continuing  effort  at  MIT  to  validate 
computer  analysis  of  welding.  A  summary  of  past  experimental 
studies  of  thermal  stresses  and  metal  movement  during  welding 
was  presented  by  Masubuchi^  and  is  reproduced  and  updated 
here  as  Table  I-A. 
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CHAPTER  II 

MATERIAL  CONSIDERATIONS 

II. A  GENERAL: 

The  Office  of  Naval  Research  is  sponsoring  an  ongoing 
program  to  investigate  thermal  stresses  and  distortion  in 
high  strength  steel  (HY-130)  structures.  HY-130  steel  is  a 
high  quality,  low  carbon,  quenched  and  tempered  steel  with 
a  yield  strength  of  130-150  ksi.  This  material  was  produced 
by  the  United  States  Steel  Company,  primarily  for  use  by  the 
U . S .  Navy . 

II. B  COMPOSITION: 

Acceptable  ranges  of  chemical  composition  for  HY-130 

steel  are  shown  in  Table  II-A.  These  ranges,  reported  by 
20 

Shackleton,  are  compositional  limits  of  HY-130  steel  as 
received  for  production. 

Mechanical  properties  of  HY-130  steel  in  the  "as  received" 

2 1 

quenched  and  tempered  condition,  as  reported  by  Schrodt,  are: 

Yield  Strength:  130-145  ksi 

Elongation  in  2  inches:  15% 

Reduction  of  Area:  50%  minimum  transverse 

70%  through  thickness 

Impact  requirements:  60  ft-lbf  at  0°  F 

A  Rockwell  hardness  test  conducted  on  the  material  to 

be  used  in  this  experimental  study  yielded  a  Rockwell  "C" 

Scale  Hardness  of  32. 


II. C  TEMPERATURE  DEPENDENCE  OF  MATERIAL  PROPERTIES: 

Subroutines  within  the  axi symmetric  finite  element 
program  developed  by  Muraki  for  analysis  of  thermal  stresses 
and  metal  movement  during  welding  require  that  material  pro¬ 
perties  be  specified  as  a  function  of  temperature.  A  search 
of  the  literature  failed  to  produce  detailed  temperature 
relationships  among  the  various  material  properties  required 
by  the  Muraki  program  for  HY-130  at  elevated  temperatures. 

Graphical  presentation  of  these  relationships  reported  by 
21 

Schrodt  are  reproduced  in  Figures  II-l  through  11-10.  A 
tabular  summary  of  this  information  is  included  as  Tables  II-B 
and  II-C . 

The  Muraki  program  which  is  to  be  verified  by  experi¬ 
mental  data  assumes  a  linear  relationship  between  various 
material  properties  and  temperature.  Assumed  linear  relation¬ 
ships  which  served  as  input  data  to  the  Muraki  program  are 
indicated  by  a  dashed  line  on  the  appropriate  figures. 
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FIGURE  II-l  TYPICAL  STRESS-STRAIN  DIAGRAM  FOR 
HY-130  AT  ROOM  TEMPERATURE 


FIGURE  II-2  ESTIMATED  EFFECT  OF  TEMPERATURE  ON  0.2% 
OFFSET  YIELD  STRESS  FOR  HY-130 
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FIGURE  II- 3  ESTIMATED  EFFECT  OF  TEMPERATURE  ON 
YOUNG'S  MODULUS  FOR  HY-130 


FIGURE  II- 4  TANGENT  MODULUS  CURVES  FOR  HY-130 
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FIGURE  II-5  ESTIMATED  EFFECT  OF  TEMPERATURE  ON 
TANGENT  MODULUS  FOR  HY-130 
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FIGURE  I 1-6  ESTIMATED  EFFECT  OF  TEMPERATURE  ON 
POISSON'S  RATIO  FOR  HY-130 
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FIGURE  II-7  ESTIMATED  COEFFICIENT  OF  THERMAL 
EXPANSION  FOR  HY-130 
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FIGURE  II-9  ESTIMATED  EFFECT  OF  TEMPERATURE  ON 
THERMAL  CONDUCTIVITY  OF  HY-130 
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FIGURE  11-10 


ESTIMATED  EFFECT  OF  TEMPERATURE  ON 
SPECIFIC  HEAT  OF  HY-130 
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SUMMARY  OF  MECHANICAL  PROPERTIES 
FOR  HY-130  STEEL 


TEMPERATURE 

(°F) 

a 

E 

H 

V 

68 

143 

30 

22.2 

.30 

200 

138 

29.6 

21.8 

.308 

400 

131.5 

28.8 

21.2 

.321 

600 

122 

27.7 

20.5 

.334 

800 

110 

26 

19.2 

.346 

1000 

88 

22.7 

17.0 

.359 

1200 

50 

14.0 

9.2 

.372 

1400 

25.5 

5.0 

4.9 

'  .384 

1600 

11.2 

4.0 

WWBM 

.397 

1800 

wwwwaw 

WBM 

2.2 

.410 

2000 

4.0 

2.0 

1.5 

.422 

_ 

2200 

2.0 

1.2 

.8 

■a 

2400 

1.0 

.4 

.4 

.447 

2500 

0 

0 

C 

.454 

o  =  Yield  Stress  2%  Offset  (ksi) 

E  =  Elastic  Modulus  (ksi  x  10 

H  =»  Tangent  Modulus  (ksi  x  10 

v  =  Poisson's  Ratio 


TABLE  II-3 
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SUMMARY  OF  PHYSICAL  PROPERTIES 


FOR  HY-130  STEEL 


TEMPERATURE 

(OF) 


68 

0.284 

200 

0.283 

400 

0.282 

600 

0.280 

800 

0.279 

1000 

0.278 

1200 

0.276 

1400 

0.275 

1600 

0.275 

1660 

0.275 

1800 

0.275 

2000 

0.275 

2200 

0.275 

8.0 


8.3 


8.18 


7.05 


6.92 


7.0 


.2 


.3 


.4 


19.4 


20.4 


21.8 


21.2 


19.5 


17.1 


15.0 


13.4 


13.1 


14.0 


0, 

.107 

0. 

.11 

0. 

.118 

0. 

.126 

0. 

.136 

0. 

.148 

0. 

.168 

0, 

.197 

0. 

.255 

0, 

.293 

0. 

.18 

0.158 


0.158 


0.18 


TABLE  II-C 


CHAPTER  III 


DESCRIPTION  OF  EXPERIMENTAL  PROCEDURE 
III. A  FABRICATION  OF  THE  WELD  SPECIMEN: 

A  circular  cylinder  of  high  strength  steel  was  formed 
in  order  to  test  the  MIT  computer  program  for  the  analysis  of 
heat  flow  in  a  weldment,  with  subsequent  use  of  the  tempera¬ 
ture  distributions  from  the  heat  flew  program  to  validate  the 
Muraki  axisymmetric  finite  element  program  for  analysis  of 
thermal  stress  and  metal  movement  during  welding.  The  HY-130 
circular  cylinder  was  formed  to  dimensions  suggested  in  the 
June  1977  proposal  on  Study  of  Residual  Stresses  and  Distor¬ 
tion  in  Structural  Weldments  in  High  Strength  Steels.  This 
research  proposal  was  made  from  Massachusetts  Institute  of 
Technology  to  the  Office  of  Naval  Research. 

Material  for  forming  the  cylinder  was  obtained  from  the 
Naval  Ship  Research  and  Development  Center,  Annapolis  Labor¬ 
atory,  Annapolis,  Maryland.  The  material  was  received  at 
Massachusetts  Institute  of  Technology  Welding  Laboratory  in 
the  form  of  HY-130  plates  60"  x  63"  x  3/4".  A  request  was 
made  to  Portsmouth  Naval  Shipyard  in  Portsmouth,  New  Hampshire, 
to  roll  the  HY-130  plates  into  cylindrical  shapes;  Portsmouth 
Naval  Shipyard  being  the  only  industrial  activity  in  the 
vicinity  of  MIT  with  experience  in  forming  HY-130.  Detailed 
drawings  and  specific  requirements  were  submitted  to  Portsmouth 
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Naval  Shipyard  and  an  agreement  was  reached  for  the  shipyard 
to  form  the  plates  into  cylinders.  However,  due  to  the 
thickness  (3/4")  of  the  plate  to  be  formed  to  a  small  outside 
diameter  (18") ,  Portsmouth  Naval  Shipyard  engineers  suggested 
that  it  would  be  necessary  to  form  a  cylinder  by  bending  two 
plates  to  semi-circular  cylindrical  shapes  and  then  welding 
the  two  halves  together.  In  view  of  this  requirement  and  in 
order  to  facilitate  easier  handling  of  the  HY-130  plates,  the 
two  plates  were  flame  cut  to  four  plates  60"  x  31V'  x  3/4" 
at  the  MIT  Welding  Laboratory  prior  to  transporting  to 
Portsmouth,  New  Hampshire. 

Stages  of  cylinder  fabrication  are  depicted  in  Figures 

III-l  and  III-2.  In  order  to  prevent  heat  resulting  from 

the  flame  cutting  operation  from  affecting  thermal  properties 

22 

of  the  steel  close  to  the  proposed  weld,  Tsai  suggested  that 
flats  required  for  bending  be  removed  by  flame  cutting  no 
closer  than  V  to  the  proposed  leading  edges  of  the  semi¬ 
circular  cylindrical  mating  surfaces.  The  flats  were  removed 
in  this  manner  and  remaining  material  to  be  removed  in  order 
to  form  a  60°  Vee  groove  with  1/8"  root,  as  illustrated  in 
Figure  III-2,  was  removed  by  mechanical  grinding.  The  semi¬ 
circular  cylinders  were  then  transported  back  to  the  MIT 
Welding  Laboratory  for  experimental  welding. 
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STAGES  OF  CYLINDER  FABRICATION 


FIGURE  III-l 


Stags  3 


Stage  h 


Stage  5 


-  Plate  halves  polled  to  semi¬ 
circular  cylindrical  shapes 
■with  flats  approximately  1  3/U' 
remaining  on  leading  edges. 


Flats  removed  by  flame  cutting 
no  closer  than  1/U”  to  mating 
edges  of  semi-cylinders. 


l/li"  excess  material  removed 
from  mating  edges  of  semi¬ 
cylinders  and  30°  angle  pro¬ 
duced  by  mechanical  grinding. 


FIGURE  III-2 
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III.B  SELECTION  OF  WELDING  CONSUMABLES: 

The  type  of  welding  suggested  in  the  research  proposal 
for  girth  welding  of  cylindrical  shells  was  multi-pass  gas 
metal  arc  welding  CGMAW) .  This  type  of  welding  requires  the 
use  of  shielding  gas  and  filler  wire,  both  of  which  can 
significantly  affect  weld  quality.  The  most  common  type  of 
shielding  gas  used  in  GMAW  of  low  carbon  steel  is  Argon  with 
25%  C02.  However,  in  the  case  of  HY-130,  experience  has 
shown  tnai  the  use  of  Argon  with  25%  C02  leads  to  inadequate 
fracture  toughness  properties  in  the  weld  deposit.  Shackleton 
recommends  rich  Argon  with  less  than  5%  Oxygen  or  10%  CC>2  to 
avoid  deterioration  in  the  fracture  toughness  properties  of 
the  weld  deposit  caused  by  Oxygen  pickup  from  the  shielding 
gas.  Shielding  gas  selected  and  purchased  for  this  experi¬ 
ment  was  Argon  with  1%  02 .  This  gas  was  reported  by 
23 

Radziminski  as  the  gas  most  effective  in  shielding  HY-130 
welds . 

Filler  wire  selected  for  the  HY-130  cylindrical  seam 
and  girth  welds  was  1/16"  diameter  Lir.de  140S  wire.  This 
wire  was  developed  by  Union  Carbide  specifically  for  HY-130 
weldments  and  was  unanimously  endorsed  in  the  literature.  Dur 
ing  periods  when  the  filler  wire  was  not  in  use,  it  was  stored 
in  a  drying  oven  containing  desiccative  material.  This  step 
was  viewed  as  necessary  in  order  to  prevent  moisture  or  other 
impurities,  or  oxides,  from  collecting  on  the  wire  surface. 


IL- 
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III.C  SELECTION  OF  WELDING  PARAMETERS : 

1.  The  variations  in  weldment  geometry  which  could 

be  considered  for  this  experiment  were  somewhat  restricted 

by  the  fact  that  cylinder  size  required  that  the  weld  be 

made  only  from  the  outside.  Based  upon  several  previous 

2  4 

welding  experiments  reported  by  Doty  and  the  above  men¬ 
tioned  restriction,  a  single  60°  Vee  groove  was  selected  as 
near  optimum  for  producing  the  seam  and  girth  welds  of  the 
cylinder.  Weld  joint  geometry  for  this  experiment  is 
depicted  below  in  Figure  III-3. 


FIGURE  III-3  WELD  JOINT  GEOMETRY  FOR  PRODUCING 
SEAM  WELDS  ON  HY-130  CYLINDER 
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2.  For  the  purposes  of  this  study,  heat  input  (Q)  is 
defined  in  the  conventional  way  as  the  product  of  arc  voltage 
and  welding  current  divided  by  arc  travel  speed,  as  expressed 
below  in  Equation  3-1 : 

HEAT  INPUT  =  Q  =  ' 1  (Equation  3-1) 

There  were  some  variations  in  welding  literature  regard¬ 
ing  the  upper  level  of  heat  input  which  would  produce  a  good 

20 

quality  weld  in  3/4"  HY-130  steel.  Shackleton  recommended 

a  maximum  heat  input  of  40,000  Joules/inch.  This  upper  limit 

25  26  23 

was  also  recommended  by  Nunez  and  Conner!  Radziminski 

recommended  a  maximum  of  45,000  Joules/inch.  Although  there 

0 

was  some  variation  in  the  suggested  upper  limits  of  heat 
input,  there  seemed  to  be  general  agreement  that  40,000  Joules/ 
inch  was  a  reasonable  value  for  the  production  of  good  quality 
welds  and  this  value  was  selected  for  production  of  the  seam 
and  girth  welds  on  the  HY-130  cylinder.  Maximum  and  minimum 
allowable  values  of  pre-heat  and  interpass  temperature  were 
also  discussed  in  the  above  quoted  references;  general  recom¬ 
mendations  were  more  varied,  but  based  on  a  rough  average  of 
recommended  values,  a  common  pre-heat  and  interpass  tempera¬ 
ture  of  200  °F  was  selected. 

3.  Filler  wire  feed  rate  is  automatically  controlled 
by  the  GMAW  welding  machine.  A  wire  feed  rate  of  approxi¬ 
mately  07  inches/minute  of  the  1/16"  Linde  140S  filler  wire 
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was  measured  during  the  production  of  a  seam  weld.  The  flow 
meter  for  the  Argon/1%  0 2  shielding  gas  was  set  at  40  ft3/hr. 
A  summary  of  welding  parameters  which  remained  unchanged  dur¬ 
ing  the  course  of  this  experiment  is  shown  in  Table  III-A 
below.  The  only  parameters  which  were  varied  between  welding 
passes  were  the  arc  voltage,  welding  current,  and  in  one 
instance  the  arc  travel  speed.  Voltage,  current,  and  arc 
travel  speed  were  recorded  for  each  pass  and  are  included  in 
the  presentation  of  results . 


WELDING  PARAMETERS 

JOINT  GEOMETRY  . 

60°  Single  V 
Groove 

MAXIMUM  HEAT  INPUT  . 

40,000  Joules/in 

FILLER  WIRE  . 

1/16"  Diameter 
Linde  140S 

feed  rate  - 

87  in/min 

SHIELDING  GAS  .  Argon/1 %  0- 

2 

flow  rate  -  40  ft3/hr 


TABLE  III-A 


SHIELDING  GAS 
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III . D  ARRANGEMENT  OF  APPARATUS: 


Welding  equipment  used  for  the  entire  experiment 
consisted  of  a  Linde  SVI-300  Power  Supply  with  ST-2  Torch 
and  associated  arc  travel  and  wire  feed  mechanisms.  Seam 
welds  were  produced  by  placing  the  cylinder  on  an  adjustable 
welding  table  and  moving  the  welding  torch  along  the  groove 
of  the  seam  weld.  Prior  to  positioning  the  cylinder  on  the 
welding  table,  the  two  cylinder  halves  were  tack  welded  at 
positions  6”  from  either  end  and  at  the  center  with  1"  welds 
on  the  inside  of  each  mating  edge.  Positions  of  the  cylinder 
and  welding  torch  for  production  of  a  seam  weld  are  depicted 
in  the  photograph.  Figure  III-4,  below. 
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III.E  TESTING  OF  WELDING  EQUIPMENT; 


In  order  to  determine  approximate  temperature  distri¬ 
butions  during  actual  welding  and  to  facilitate  a  thorough 
checkout  of  welding  equipment  prior  to  welding  passes  for 
data  gathering,  it  was  deemed  advantageous  to  produce  a  seam 
weld  on  one  side  of  the  cylinder  without  benefit  of  instru¬ 
mentation  to  record  temperature  and  strain.  This  allowed 
the  wire  feed  mechanism,  arc  travel  controls  and  indicators, 
gas  flow  regulators  and  meters,  and  the  welding  equipment 
electrical  circuity  and  instrumentation  to  be  tested  and 
checked  for  accuracy  under  conditions  expected  during  the 
actual  test.  This  "test  run"  also  provided  an  opportunity 
for  operator  training.  A  preliminary  checkout  of  welding 
equipment  had  been  conducted  using  "bead-on-plate"  welds  on 
scrap  pieces  of  HY-130. 

Approximate  temperature  distributions  to  be  expected 
during  the  instrumented  test  were  determined  by  the  applica¬ 
tion  of  temperature  indicating  liquid  (TEMPILAQ)  in  narrow 
bands  transverse  to  the  welding  direction.  TEMPILAQ  liquids 
corresponding  to  temperature  levels  of  400  °F,  600  °F,  800  °F, 
1000  °F,  and  1100  °F  were  applied  to  the  chemically  cleaned 
inner  and  outer  surfaces  of  the  cylinder.  TEMPILAQ  was  re¬ 
applied  between  passes  and  distances  of  temperature  penetra¬ 
tion  from  the  weld  centerline  on  both  the  inside  and  outside 


i 


45 


of  the  cylinder  were  recorded  for  each  welding  pass.  Welding 
current,  arc  voltage,  arc  travel  speed,  and  maximum  distances 
from  the  weld  centerline  for  the  various  TEMPILAQ  levels  were 
recorded  for  each  welding  pass  and  are  shown  in  Table  III-B. 
This  tabulated  data  was  used  as  a  basis  for  determining  strain 
gage  and  thermocouple  locations  for  the  instrumented  seam  and 
girth  welds  to  be  used  in  validating  heat  flow  and  thermal 
stress  computer  programs. 
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WEEDING  PASS 

♦TEMPERATURE  PENETRATION  DISTANCE 

1.  I  «  26QA 

7  =  23V 
t  —  12  in/min 

OUTSIDE  INSIDE 

LEFT  RIGHT  LEFT  RIGHT 

UOOF  .86  -  ,9U  .88 

600F  -  -  .63  .60 

8 OOF  -  .UU  .ho 

1000F  -  .38  .35 

HOOF  -  -  .25  .22 

2.  I  a  260A 

7  *  23V 
▼  *  12  in/min 

UOOF  .61  .38  .6 9  .75 

600?  -  -  .38  .56 

8 OOF  -  .22  .38 

1000F  -  .06  .31 

HOOF  -  .19 

3.  1  *  260A 

V  -  237 
▼  a  12  in/min 

UOOF  .80  -  .88  .63 

600F  -  -  .50  .38 

800?  -  .13 

10Q0F  -  - 

HOOF  -  - 

U.  Ia  26QA 

7  a  23V 
▼  *  12  in/min 

UOOF  -  .92  .50  .88 

600F  -  .U9  .19  .50 

8 OOF  -  - 

1000F  -  - 

HOOF  -  - 

5.  I  a  35CA 

7  =  227 
v  a  10  in/min 

UOOF  .92  1.55  1.06  1.63 

600F  -  .80  .38  .75 

800F  -  .58  -  .38 

1000F  .09  - 

HOOF  .03  - 

6.  I  ■  33QA 

7  =  217 
v  -  10  in/  min 

UOOF  1.U9  1.15  1.50  .  88  ! 

600F  .80  .75  .75 

8qof  .55  .25  i 

looo?  ,U9  j 

HOOF  .U9  ; 

7.  I  =  35QA 

V  a  227 
▼  •  20  in/min 

1 

UOOF  .U6  .78  .38  .75 

6oof  -  .UO  -  .13 

800?  .25  - 

100  OF  .03  - 

HOOF  -  - 

8.  I  =  370A 

7  =  217 
v  *  20  in/min 

TEMPERATURE  PENETRATION  DISTANCES  NOT  RECORDED 

*  Temperature  penetration  distances  are  expressed  in  inches,  measured 
in  the  direction  indicated  from  the  weld  centerline. 


TABLE  III-B 


III.F  SELECTION  AND  ARRANGEMENT  OF  INSTRUMENTATION  TO 

MEASURE  STRESS  AND  DISTORTION  DURING  WELDING: 

Stresses  occurring  during  the  thermal  cycles  associated 
with  seam  and  girth  welding  of  the  HY-130  cylinder  were  de¬ 
termined  by  measuring  strains  as  a  function  of  time  (i.e.  - 
location  of  the  strain  measurement  with  respect  to  the  time 
history  of  welding  arc  location)  and  as  a  function  of  distance 
from  the  weld  centerline.  Strains  were  measured  by  placing 
electrical-resistance  wire  gages  on  the  outside  surface  of 
the  cylinder.  Bi-axial  strain  gages  were  placed  with  one 
axis  aligned  with  the  cylinder  axis  and  welding  direction 
(herein  called  the  longitudinal  direction)  and  the  other  axis 
aligned  in  a  direction  normal  to  the  welding  direction  (herein 
referred  to  as  the  transverse  direction.  Characteristics  of 
the  strain  gages  used  for  this  experiment  are  shown  in 
Table  III-C. 

In  order  to  validate  the  MIT  heat  flow  program,  it  was 
necessary  to  measure  the  heat  distribution  within  the  cylinder. 
Heat  flow  and  temperature  distribution  within  the  HY-130 
cylinder  was  viewed  as  critical  to  validation  of  the  Muraki 
program,  inasmuch  as  strain,  stress,  and  distortion  calcula¬ 
tions  are  all  based  upon  temperature  distribution  data  which 
must  be  supplied  as  input  to  the  program.  Thermocouples  used 
for  measuring  temperature  distribution  were  Chromel/Alumel 
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STRAIN  GAGE  CHARACTERISTICS 


Type  Used:  FAET-12D-12S6 (EL) 

Manufacturer:  BLH  Electronics,  Inc. 


Waltham,  Massachusetts  02154 


TABLE  III-C 


Strain  indications  were  corrected  for  temperature 
by  the  following  temperature  error  (apparent  strain)  equation 

EAP  =  -91,37  +  2.17T  -  2.33  x  1Q-2T2  +  5.31  x  10~5T 

-a  4 

-  2.40  x  10  T 
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wires  which  were  peered  into  small  (0.040")  indentations 
drilled  into  the  cylinder  outer  surface.  Thermocouple 
circuit  continuity  was  assured  by  the  use  of  an  inductive 
spot  welder  on  the  thermocouple  connections.  After  circuit 
continuity  was  assured,  the  thermocouples  were  protected  by 
the  use  of  sealing  cement. 

Since  longitudinal  and  transverse  strains  were  to  be 
measured  with  gages  which  required  temperature  correction, 
it  was  necessary  to  place  thermocouples  in  positions  which 
would  give  temperature-time  histories  consistent  with  strain 
gage  locations.  Arrangement  of  strain  gages  and  thermo¬ 
couples  on  the  HY-130  cylinder  outer  surface  is  shown  in 
Figure  III-5.  Arranging  the  strain  gages  and  thermocouples 
in  this  manner  allowed  them  to  be  used  for  measuring  strain 
and  temperature  distributions  during  the  production  of  the 
second  seam  weld  (the  first  seam  weld  was  made  without  in¬ 
strumentation)  and  then  to  be  used  again  after  the  cylinder 
had  been  cut  into  sections  and  mating  edges  prepared  for  a 
girth  weld.  Arranging  the  strain  gages  and  thermocouples 
at  a  45°  angle  with  respect  to  the  seam  weld  and  at  equi¬ 
distant  locations  with  respect  to  each  other  allowed  the 
strain  gages  and  thermocouples  to  give  strain-temperature 
distributions  which  could  be  adjusted  in  time  so  as  to  be 
coincident . 
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It  should  be  pointed  out  that  the  practice  of  locating 
strain  gages  and  thermocouples  so  as  to  require  time  adjust¬ 
ment  of  their  readings  assumes  quasi-steady-state  heat  flow 
(i.e.  -  the  heat  flow/temperature  distribution  with  respect 
to  the  welding  arc  is  the  same  at  a  position  40"  along  the 
length  of  the  cylinder  as  at  a  position  41"  along  the  cylin¬ 
der)  .  This  assumption  is  necessary  in  view  of  the  fact  that 
strain  gages  and  thermocouples  cannot  be  located  at  the  same 
position . 

Arranging  the  strain  gages  and  thermocouples  at  a 
45°  angle  with  respect  to  the  seam  weld  was  viewed  as  ad¬ 
vantageous  since  this  scheme  places  the  strain  gages  and 
thermocouples  which  were  subjected  to  the  most  severe  thermal 
cycle  during  production  of  the  seam  weld  in  a  position  which 
is  subjected  to  the  least  severe  thermal  cycle  during  the 
girth  weld.  Failures  occurring  during  the  seam  weld  are 
most  likely  to  occur  at  locations  which  are  least  critical 
to  measurements  to  be  taken  during  the  girth  weld. 

After  determining  and  marking  the  desired  strain  gage 
locations,  strain  gages  were  cemented  to  the  cylinder  outer 
surface.  Physical  size  and  weight  of  the  cylinder  prevented 
the  usual  technique  of  curing  the  strain  gage  bonding  cement 
in  laboratory  ovens.  In  order  to  assure  that  the  strain  gage 
bonding  would  withstand  expected  temperatures,  a  5  foot  by 
5  foot  thermal  insulating  blanket  (graciously  manufactured 
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for  this  purpose  by  Bethlehem  Steel  Corporation  Shipyard  in 
Boston,  Massachusetts)  was  placed  around  the  cylinder  and 
electrical  heating  elements  placed  inside.  The  temperature 
of  the  cylinder  was  slowly  raised  to  250  °F  and  the  strain 
gage  bonding  cement  was  cured  at  this  temperature  for  a  period 
of  four  hours,  as  directed  by  the  strain  gage  manufacturer's 
instructions.  The  heating  elements  used  for  this  purpose 
were  later  used  to  preheat  the  cylinder  prior  to  welding. 

For  production  of  the  girth  weld,  it  was  necessary  to 
cut  the  HY-13G  cylinder  into  sections.  Original  plans  called 
for  cutting  the  60  inch  cylinder  into  two  30  inch  sections 
which  would  have  mating  edges  prepared  for  a  girth  weld  by 
machining  a  60°  single  Vee  groove.  Due  to  physical  limitations 
of  laboratory  facilities  it  was  necessary  to  reduce  the  length 
of  the  two  sections  to  be  joined  by  a  girth  weld  to  20  inches 
each.  This  required  that  the  HY-130  cylinder  be  cut  into 
three  sections,  each  20  inches  in  length.  This  decision  was 
made  prior  to  production  of  the  instrumented  seam  weld,  and 
the  strain  gages  were  placed  so  as  to  permit  measurements 
during  the  seam  weld  as  well  as  during  production  of  the  girth 
weld. 

After  completion  of  the  seam  weld,  the  60  inch  long 
HY-130  cylinder  was  transported  to  the  Army  Materials  and 
Mechanics  Research  Center  in  Watertown,  Massachusetts,  where 
it  was  sawed  into  three  sections,  each  20  inches  in  length. 


A  cross-section  of  the  previously  produced  seam  weld 
was  removed  from  the  h  inch  thick  ring  and  later  ground, 
polished,  and  etched  for  analysis  of  the  metallurgical  proper 
ties  of  the  seam  weld.  Details  of  this  procedure  and  a  dis¬ 
cussion  of  the  results  are  contained  in  Chapter  IV  of  this 
study . 

Mating  edges  of  the  two  sections  to  be  joined  by  girth 
welding  were  then  prepared  for  welding  by  machining  a  60° 
single  Vee  groove  as  shown  below  in  Figure  III-7.  Heat 
buildup  during  machining  was  prevented  by  the  use  of  a  con¬ 
tinuously  flowing  oil  and  water  cooling  bath  during  machining 


included 


3/4  inch 


FIGURE  III-7 
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The  two  20  inch  cylindrical  sections  were  transported 
back  to  the  MIT  welding  laboratory  where  apparatus  required 
for  production  of  the  girth  weld  was  assembled.  In  order  to 
control  arc  travel  speed  during  production  of  the  girth  weld, 
it  was  necessary  to  design  and  manufacture  a  supporting  struc¬ 
ture  for  the  cylinder  which  would  permit  cylinder  rotation  at 
a  controlled  speed  with  respect  to  the  fixed  position  of  the 
GMAW  torch.  This  was  accomplished  by  producing  a  structure 
which  was  designed  to  support  the  two  20  inch  cylinders  from 
within.  Details  of  this  supporting  structure  are  shown  below 
in  Figure  III-8. 


FIGURE  III-8 
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The  rectangular  base  of  the  supporting  structure  was 
bolted  to  the  base  of  a  tilting  turntable  and  the  20  inch 
long  cylinder  halves  were  fixed  to  the  support  by  force 
exerted  on  the  interior  walls  of  the  cylinder  by  the  threaded 
dowels.  The  turntable  was  then  tilted  so  as  to  place  the 
now  40  inch  long  cylinder  (with  60°  Vee  groove  located  in  the 
center)  in  a  horizontal  position.  The  threaded  dowels  were 
adjusted  in  order  to  provide  for  constant  clearance  between 
the  stationary  GMAW  torch  and  cylinder  outer  wail  during 
cylinder  rotation.  Turntable  controls  were  adjusted  and 
calibrated  in  order  to  achieve  desired  arc  travel  speeds. 

A  photograph  of  the  apparatus  in  use  during  production  of  the 
girth  weld  appears  below  as  Figure  III-9. 


FIGURE  III-9 
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In  order  to  measure  radial  distortion  during  production 
of  the  girth  weld,  two  Bourns  Linear  Motion  Potentiometers 
were  placed  at  the  inner  surface  of  the  cylinder.  Character¬ 
istics  of  the  linear  motion  potentiometers  are  shown  below 
in  Table  III-D. 


ACCEPTANCE  TEST  RECORD 
FOR  LINEAR  MOTION  POTENTIOMETERS 


TESTS  PERFORMED  SPEC  LIMITS 


MECHANICAL  TRAVEL  . 

,  . .  0 . 44  in  max. 

. 

i 

ACTUATION  FORCE  . . 

.  .  .  4  oz  min 

TOTAL  ELEMENT  RESISTANCE  . . . 

...  4750  to  5250  Ohms 

J 

1 

INSULATION  RESISTANCE  _ 

.  .  .  50  Megohms  min 

at  500  V  DC 

j 

NOISE  OR  CONTACT 

j 

RESISTANCE  VARIATION . 

100  ohms  max 

1 

END  SETTINGS  . 

, .  0 . 1  to  5 . 0% 

each  end 

i 

LINEARITY  . 

.  *  1.0%  VR  max 

TABLE  III-D 
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Displacement  transducer  circuitry  and  arrangement  cf 
linear  motion  potentiometers  within  the  HY-130  cylinder  are 
indicated  in  Figures  III-10  and  III-ll  respectively.  Prior 
to  mounting  the  displacement  transducers  in  position  within 
the  cylinder,  a  calibration  test  was  conducted  by  deflecting 
the  displacement  transducers  through  measured  distances  and 
recording  resultant  voltages.  A  plot  of  this  calibration 
test  appears  as  Figure  III-12.  After  completion  of  calibra¬ 
tion  tests,  the  linear  motion  potentiometers  were  installed 
within  the  cylinder  as  indicated  in  Figure  III-ll.  Holes 
were  drilled  and  threaded  in  the  cylinder  support  mechanism 
and  the  linear  motion  potentiometer  support  reds  were 
threaded  and  screwed  into  these  holes  at  diametrically 
opposed  locations  within  the  cylinder.  Potentiometers 
corresponding  to  Channels  1  and  2  of  the  recorder  were 
located  at  positions  1/2  inch  and  1  inch  respectively  in 
the  cylinder  axial  direction  from  the  weld  centerline. 
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CHANNEL  1  DISPLACEMENT  TRANSDUCER 
SENSING  ELEMENT  LOCATED  H  INCH 
FROM  WELD  CENTERLINE.  INSULATED 
FROM  CYLINDER  INNER  SURFACE 
BY  CERAMIC  BEAD. 


DISPLACEMENT  TRANSDUCER  SENSING  E 
MOUNTING  RODS,  SCREWED  INTO 
CYLINDER  SUPPORT. 


CHANNEL  2  DISPLACEMENT  TRANSDUCER 
SENSING  ELEMENT  LOCATED  1  INCH  FROM 
WELD  CENTERLINE.  ALSO  INSULATED. 


CYLINDER  SUPPORT  MECHANISM 


FIGURE  III-ll 


DISPLACEMENT  TRANSDUCER  SENSING  ELEMENT 
LOCATIONS 
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Strain  and  temperature  were  measured  at  the  indicated 
locations  during  production  of  the  instrumented  seam  weld; 
and  strain,  temperature,  and  distortion  were  measured  during 
production  of  the  girth  weld.  Equipment  used  to  record  the 
measured  values  included  the  following: 

Minneapolis-Honeywell  H1508  2  4  Channel  Visicorder 

Three  Helland  Universal  Bridge  3alancer  Units, 

Model  no.  82-6  (ser  no's  1332,1307,1312) 

16  Channel  Strain  Gage  Calibrator 
(one  arm  bridge  adapter) 

Gould  Brush  220  Two-Channel  Recorder  (distortion) 

General  arrangement  and  operation  of  this  equipment 
during  production  of  the  girth  weld  are  indicated  below  in 
Figure  III-13. 


FIGURE  III-13 
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III.G  PRODUCTION  OF  THE  INSTRUMENTED  SEAM  AND  GIRTH  WELDS : 

1.  The  Seam  Weld:  After  the  HY-130  cylinder  had  been 
suitably  supported  and  aligned;  instrumentation  and  other 
apparatus  arranged,  checked  and  calibrated,  the  instrumented 
seam  weld  was  produced  by  making  seven  passes.  Figure  III-14 
below  shows  approximate  dimensions  and  location  of  each  pass. 
Values  of  weld  parameters  for  each  pass  are  tabulated  in 
Table  III-E. 


( 


FIGURE  III-14 
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2.  The  Girth  Weld:  Production  of  the  girth  weld 
required  careful  positioning  of  the  rotating  HY-130  cylinder 
with  respect  to  the  GMA  welding  torch.  The  girth  weld  required 
six  passes  with  approximate  dimensions  shown  below  in  Figure 
III-15.  Girth  weld  parameters  are  also  included  in  Table 


III-E. 


I  =  Welding  Current  in  Amperes 
V  =  arc  voltage  in  Volts 

4 

v  =  arc  travel  speed  in  inches/sec 


TABLE  1'II-E  SEAM  AND  GIRTH  WELD  PASSES 
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CHAPTER  IV 

METALLURGICAL  CONSIDERATIONS 
IV. A  PREPARATION  OF  SPECIMEN: 

When  the  60  inch  cylinder  was  sawed  into  sections,  as 
previously  discussed  in  Chapter  III,  a  cross-section  of  the 
instrumented  weld  was  cut  from  the  %  inch  ring.  The  weld 
specimen  was  cut  to  approximately  3/4  inch  by  3/4  inch  and 
marked  for  identification  as  shown  in  Figure  IV- 1.  Mounting 
for  polishing  was  in  thermo-setting  bakelite  resin,  pressed 
into  a  1**  inch  diameter  by  3/4  inch  long  cylindrical  shape 
under  hydraulic  pressure  of  4200  psi  at  250  °F  for  2  minutes 
using  a  Buehler  press. 

The  specimen  was  then  rough  polished  by  successive  use 
of  60,  120,  180,  240,  320,  and  400  grit  3M  TRI-M-ITE  Silicon 
Carbide  paper.  Fine  polishing  was  accomplished  by  successive 
use  of  l.Qu,  0.3y,  and  0.0 5y  alumina  polishing  powder.  Etching 
for  photo-micrographs  required  about  15  seconds  with  5%  NITAL. 
The  5Q0X  photo-micrographs,  Figures  IV -2  through  IV-5 ,  were 
produced  by  the  use  of  a  LEITZ  WETZLAR  microscope  with  NPL 
50X/0.85  lens,  using  Polaroid  type  52  POLAPAN  film.  Fine 
polishing  was  repeated  and  the  specimen  was  etched  using  1% 
NITAL  prior  to  the  production  of  1200X  and  6000X  electron  scan 
micrographs.  These  electron  scan  micrographs  appear  as  Figures 
IV- 6  through  IV-13.  Locations  of  the  photomicrographs  and  the 
electron  scan  micrographs  are  indicated  on  Figure  IV-1. 


FIGURE  IV- 1  DIAGRAM  SHOWING  LOCATIONS  OF  PHOTO-MICROGRAPHS 
AND  ELECTRON  SCAN  MICROGRAPHS.  3/4"  SAMPLE 
WAS  CUT  FROM  THE  SEAM  WELD  (SHOWING  CROSS 
SECTION  OF  WELD  -  WITH  NOTCH  IN  UPPER  LEFT 
CORRESPONDING  TO  WELD  CENTERLINE) .  AFTER 
POLISHING  AND  ETCHING, 50 OX  PHOTO-MICROGRAPHS 
AND  120 OX/6 00 OX  ELECTRON  SCAN  MICROGRAPHS  WERE 
PRODUCED  AT  THE  LOCATIONS  INDICATED. 
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FIG  IV-4  PHOTOMICROGRAPH  OF  HY-130 
HAZ.  500X .  5%  NITAL. 


FIG  IV- 5  PHOTOMICROGRAPH  OF  HY-130 
PARENT  METAL.  500X.  5%  NI 


G  IV- 6  ELECTRON  SCAN  MICROGRAPH  OF 
HY-130  WELD  METAL.  1200X. 
1%  NITAL. 


FIG  IV-7  ELECTRON  SCAN  MICROGRAPH  OF 
HY-130  WELD  METAL.  6000X. 
1%  NITAL. 
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FIG  IV- 10  ELECTRON  SCAN  MICROGRAPH  OF 

KY-130  HA2 .  1200X .  1%  NITAI 


j#$ 
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FIG  IV- 11  ELECTRON  SCAN  MICROGRAPH  OF 
HY-130  HAZ.  6000X.  1%  NITAL 


FIG  IV- 1 2  ELECTRON  SCAN  MICROGRAPH  OF 
KY-I30  PARENT  METAL.  1200X. 
1%  NITAL. 
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IV- 13  ELECTRON  SCAN  MICROGRAPH  OF 
HY-130  PARENT  METAL .  6000X 

1%  NITAL. 
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CHAPTER  V 

DESCRIPTION  OF  ANALYTICAL  TOOLS  AND  PROCEDURES 
V . A  INTRODUCTION: 

Analytical  procedures  used  in  this  study  to  predict 
temperature  distribution,  strain,  stress,  and  distortion 
occurring  during  production  of  a  weldment  were  developed  at 
Massachusetts  Institute  of  Technology.  Both  the  heat  flow 
program  and  the  program  used  to  analyze  thermal  stress  and 
metal  movement  use  finite  element  techniques  and  are  written 
in  the  FORTRAN-IV  language. 

Discussion  of  these  programs  within  this  study  will  be 
limited  to  a  brief  synopsis  of  their  theoretical  background, 
presentation  of  assumptions  and  applicable  equations  upon 
which  they  are  based,  and  a  summary  of  welding  parameters 
pertinent  to  production  of  the  HY-130  weldment  which  were 
used  as  input  to  the  programs.  Detailed  discussion  of  these 
programs,  lists,  card  decks,  and  manuals  describing  their  use 
are  located  within  the  Department  of  Ocean  Engineering  at 
Massachusetts  Institute  of  Technology. 

V.B  DESCRIPTION  OF  THE  HEAT  FLOW  PROGRAM: 

1.  Theoretical  Background.  The  finite  element  computer 

program  used  in  this  study  to  generate  temperature  distribution 

in  the  HY-130  cylinder  girth  weld  is  Program  B-4  in  a  report 

2  7 

developed  by  Muraki  and  Masubuchi  entitled  "Computer  Programs 
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Useful  for  the  Analysis  of  Heat  Flow  in  Weldments"  at  MIT  in 
1974.  This  program  develops  temperature  distribution  based 
on  heat  conduction  in  quasi-stationary  state  from  a  moving 
point  heat  source  in  a  plate  of  finite  thickness.  The 
solution  assumes: 

a.  The  size  of  the  heat  source  is  small  compared 
to  that  of  the  body  under  consideration. 

b.  Heat  is  supplied  instantaneously. 

c.  Thermophysical  properties  do  not  depend  upon 
temperature . 

d.  Heat  is  transferred  in  conduction  only  (other 
modes  such  as  radiation  and  convection  are  not  considered) . 

e.  The  material  of  the  body  is  uniform  and  iso¬ 
tropic  in  the  thermophysical  properties. 

The  solution  is  based  upon  the  following  equation  which 
expresses  the  difference  between  the  temperature  at  time  t 
(in  seconds  after  the  arc  is  initiated)  and  the  initial 
temperature : 


T-T  = 

o  2~tc 


Where : 


-VS 

-vr 

-vrA\ 

~~ 

e  2K 

e  2K 

+  /  |e  2K 

+  e  2K  1 
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\  n 
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(Eq.  5-1) 


"a 

x  - 

vt 

2 

2 

r* 

=  X 

+  y 

2 

r 

=  s 

+  y 

n 

2 

r 1 
n 

=  S 

+  y 

+  z 


+  y  +  (2nh  -  z) 
+  y2  +  (2nh  +  z) 
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2 .  Geometry  of  Weld  Specimen  and  Description  of 
Finite  Element  Representation  of  Temperature  Distribution 


Program: 


Division  of  the  HY-130  cylinder  girth  weld  into  elements 
for  use  in  the  heat  conduction  program  is  illustrated  in 
Figure  V-l.  Node  coordinates  (x,  y,  and  z  values)  for  each 
of  the  114  nodes  are  listed  in  Table  V-A.  The  weldment  is 
represented  by  finite  elements  which  are  arranged  so  as  to 
cause  node  locations  for  analytical  determination  of  heat 
distribution  to  be  coincident  with  locations  of  thermocouples 
which  were  actually  placed  on  the  cylinder  outer  surface  in 
order  to  experimentally  determine  temperature  distribution. 


Node  coordinates  listed  in  Table  V-A  are  based  upon  an 
origin  (x~y=z=0)  which  is  coincident  with  the  starting  point 
of  the  heat  source  (welding  arc) ,  with  subsequent  times  of 
interest  corresponding  to  location  of  the  heat  source  (in 
nodal  coordinates)  at  the  referenced  time.  This  analysis 
assumes  that  the  heat  source  moves  on  the  plane  z=0  along 
the  x-axis.  The  x  direction  (direction  of  welding  arc 
travel)  would  appear  on  Figure  V-l  as  normal  to  the  y-z 
plane,  positive  direction  into  the  paper.  In  this  figure 
the  y-axis  is  along  the  cylinder  length  and  the  z-axis  is 
in  the  radial  direction. 


TE  ELEMENT  REPRESENTATION  OF  WELDMENT 
GEOMETRY  FOR  HEAT  FLOW  ANALYSIS 


NODE  COORDINATES  (HEAT  FLOW  PROGRAM) 
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In  order  to  determine  temperature  distribution  at 
weldment  nodes  for  the  first  pass  of  the  HY-130  cylinder 
girth  weld,  the  following  parameters  were  used  as  input  to 
the  heat  flow  program: 

Velocity  of  welding  arc  (v)  =  0.2000  in/sec 
Heat  intensity  (q)  =  6034  Joules/ 

I  x  V  x  n 


sec  in 


Where:  q  = 


arc 


(Eq.  5-2) 


n 


I  =  welding  current 
V  =  arc  voltage 
h  =  plate  thickness 
n  .  =  welding  arc  efficiency 

cl  IT  C 

Density  (o)  =  0.282  lbm/^3 

Specific  heat  (<)  =  124.4664  Joules/^^  o? 

Thermal  conductivity  (c)  =  0.5250  Joule/sec  c_ 

It  should  be  noted  that  this  stage  of  the  computer 

analysis  does  not  provide  for  temperature  dependence  of 

material  properties.  Values  used  in  this  case  to  determine 

temperature  distribution  are  these  valid  at  400  °F .  It  is 

also  significant  to  note  that  for  this  analysis,  an  arc 

efficiency  (n  )  of  0.65  was  assumed. 

2-1- 
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Nc  modifications  were  required  to  the  analytical  or 
mathematical  techniques  required  for  generating  temperature 
distributions  to  be  used  by  the  Muraki  strain,  stress,  and 
distortion  calculating  program.  However,  since  large  amounts 
of  numerical  data  were  generated  by  the  heat  flow  program  for 
input  to  the  Muraki  program,  it  was  necessary  to  modify  the 
heat  program  output.  Since  weldment  geometry  and  finite 
element  representation  of  the  nodes  and  elements  were  the 
same  (except  for  coordinates)  for  both  programs,  the  heat  flow 
program  was  modified  so  as  to  provide  its  output  in  the  form 
of  punched  cards  suitable  for  use  as  input  to  the  Muraki 
program.  This  reduced  the  probability  of  introducing  numeri¬ 
cal  error  between  the  two  programs  as  well  as  improving  the 
ease  of  program  interaction. 
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V.C  DESCRIPTION  OF  THE  PROGRAM  TO  CALCULATE  STRAIN,  STRESS, 
AND  DISTORTION; 

1.  Theoretical  Background:  The  program  used  in  this 
study  to  analytically  generate  strain,  stress,  and  distortion 
distributions  within  the  HY-130  cylinder  girth  weld  was  de¬ 
veloped  at  MIT  by  Toyohiki  Muraki  for  thermo  elasto-plastic 
analyses  of  axisymmetric  structural  components .  The  Muraki 
program  is  written  in  FORTRAN  IV  language  and  uses  finite 
element  techniques  in  an  incremental  fashion  to  calculate 
strains,  stresses,  and  distortion. 


The  Muraki  program  has  great  flexibility  in  terms  of 

allowing  temperature  dependency  of  material  properties  to  be 

considered  as  well  as  considering  the  effect  of  melting  and 

deposition  of  metal  which  occurs  during  welding.  Weldment 

restraint  and  concentrated  loads  are  also  considered  i..  the 

analytical  solution.  In  the  "Manual  on  Axisymmetric  Finite 

Element  Program  for  Analysis  of  Thermal  Stress  ar.d  Metal 

2  3 

Movement  During  Welding",  Muraki  presents  a  brief  synopsis 
of  the  theoretical  background  upon  which  this  analysis  is 
based.  This  synopsis  is  reproduced  here  as  follows: 


Principle  of  Virtual  Work 

The  principle  of  virtual  work  for  the  quasi-sta- 
problem  is,  in  general, 


///„{  o £  }T{  a}dV  -  <5u}T{p}dV  -  ff  {  m 


v 


V 


1  -  — .  - . 
;  F  ;  c ; 


=  0  (1) 


32 


where  superscript  T  stands  for  the  transpose  of  a  vector  or 
a  matrix.  {e},  {a},  and  {u}  are  vectors  of  the  time  rates  of 

strain,  stress  and  displacement,  respectively. 

«  * 

P  and  F  denote  the  time  rates  of  body  force  over  the 
volume,  V,  and  surface  force  on  the  boundary  s  .'face ,  S^, 
respectively . 

The  boundary  condition  on  the  specified  displacement  on 
the  surface,  S_ ,  is  given  as 

{u}  =  {u}  on  S2  (2) 

It  is  noted  that  the  problem  here  is  formulated  in  the 
small  displacement  theory  based  on  the  rectangular  cartesian 
coordinate  system. 

Constitutive  Equations 

One  of  the  important  features  in  the  analysis  of  thermal 
stress  and  distortion  during  welding  is  the  temperature  de¬ 
pendency  of  material  properties  and  the  yield  criterion. 
Therefore  the  constitutive  equations  are  provided  here  in  a 
general  form. 

The  rate  of  strain  for  the  elasto-plastic  state  is  given 

by : 

u)  =  u(e)}  +  a(p>>  i3) 

where  superscripts  e  and  p  refer  to  elastic  and  plastic  com¬ 
ponents  of  the  strain,  respectively. 
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The  rate  of  plastic  strain  is  assumed  as 

{s(p)  }  =  A(g)  (4) 

where  A  is  a  proportional  constant  and  g  is  derivatives  of 
the  yield  function,  f,  with  respect  to  the  stress, 

(if  j  =  1,2,3); 

a  =  1* 

3  3cis  (5) 

Using  Hooke's  law  and  Equation  (4),  we  have 

{e}=g-|De|1{a}+A{g}  (6) 

where  E | De ]  is  the  matrix  form  of  Hooke's  law. 


Thermal  strain  and  the  temperature  dependency  of 
material  properties  are  added  to  Equation  (6)  as  follows: 


tcr\  =  £, 


E  |  De  !  _1  {c}-  Iq-  |  De  i  ~1{c} 


-1 


A  +  s  {  0  i  j  }  +  A  { g } 


(7) 


.  a 


where  e "  is  the  thermal  strain,  and  6. .  is  the  Kronecker 
symbol . 


Equation  (7)  is  rewritten  fo: 


:e  rate  of  stress  as, 


(a}=  ”|De[  (s}  +  — { c }  -  E I  Del {6..} 

-  AE  j  De  j  {g} 


(8) 


Note : 


stands 


for  the  inverse  of 


the  matrix 
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If  the  von  Mises '  yield  criterion  is  adopted,  f  in 
Equation  (5)  is  given  by 

f  =  a  -  C  (ef?*  ,  T)  (9) 

where  C  is  related  to  the  strain  hardening  and  given  as  a 
function  of  plastic  strain,  ef?^  and  temperature,  T. 

a  is  referred  to  as  the  equivalent  stress  and  given 
by 

a2  =  3/2ic'}Tia'}  (10) 

It  is  noted  that  (a"}  is  the  vector  of  deviator  stresses 
which  is  different  from  {cr}. 


When  a  state  of  stress  in  a  portion  of  a  body  is  plastic 
and  then  changes  to  a  new  state  of  stress  which  is  still 
plastic,  the  following  relation  is  obtained  from  Equation 
(9 )  : 


_  r  -  rc ‘ )  ‘  w  1 

3a  3s^) 


0 


(11) 


Equation  (11)  with  Equations  (4),  (9) ,  and  (10)  leads 

to 

{a'}T{a}  =  -2/3a(-  AH"  +  ||  T)  (12) 

where 


From  Equations  (8)  and  (12),  A  is  obtained  as  follows: 


and 
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A  =  - i - = -  x  (E  { a  "  } T  |  De  '  {  e  } 

2/3  a  H'  +  E{a'}  |Dej {g} 

+  E  {a'}T  Ca}  -  £9E{a^}T!De| {5. .} 

E  13 

+2/3  ct  ~  T)  (13) 


Substituting  A  of  Equation  into  (3)  we  have  a  general 
form  of  the  stress-strain  relationship  as 


{a}  =  (E  |  De  !  -  1/SqE  |  De  j  {g}  E{  cr "  }T  |  De  !  )  x 
({£}  - 

+  E/E  { 1 1 1  -  1/S  E | De | {g}{c"}T) { a } 

-  l/So  E | De ! { g }  2/3c  ||  T 


(14) 


where 

SQ  =  2/3  a  E'  +  E{a"}T|De| {g} 
and  | I !  is  called  the  unit  matrix. 


Equation  (14)  covers  the  stress-strain  relationship  in 
the  elastic  as  well  as  plastic  state.  This  can  be  easily 
confirmed  by  eliminating  terms  related  to  the  yield  function 
from  Equation  (14) . 


Yield  Criterion 


As  described  in  Equation  (9)  the  yield  function  used 
here  is  that  proposed  by  von  Mises. 

The  yield  function  is  to  separate  the  region  of  elastic 


deformation  from  that  of  elastic  deformation.  Since  the 


deformation  is  closely  related  to  the  loading  process,  the 
yield  function  is  also  called  a  loading  surface  which  is 
defined  in  stress  space. 

If  a  part  of  a  material  in  question  is  in  the  region  of 
plastic  deformation,  which  means  that  the  state  of  stress  of 
the  part  is  located  on  the  loading  surface,  one  of  the 
following  states  will  be  realized  at  the  next  instance: 

• 

f  =  0,  S'  <  0  (unloading) 

f  =  0,  S'  =  0  (neutral) 

* 

f  =  0,  S'  >  0  (loading) 

where 

f  =  .  if  t 

3 c ij  +  3T 

ij 

If  a  part  of  the  material  is  in  the  region  of  elastic 
deformation,  one  of  the  following  states  will  be  realized 
at  the  next  instance: 

f  <  0  (elastic) 

(16) 

f  =  0  (plastic) 

It  is  worth  looking  into  Equation  (15)  more  closely. 
Using  Equations  (10)  and  (12)  with  Equation  (15),  the 
following  relationship  is  obtained: 

S'  =  AH'  (17) 

Equation  (17)  indicates  that  when  non-strain  hardening 


is  assumed  of  a  material,  Equation  (15)  cannot  be  used  to 
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judge  whether  or  not  a  part  of  the  material  recovers  from 
the  plastic  to  the  elastic  state  of  stress.  To  avoid  this 
difficulty,  Equation  (17)  will  be  written  in  the  form  of: 


iim  _ 
H^O  H»  ~ 


(18) 


Therefore  Equation  (15)  is  written  as, 
f  =  0,  A  <  0  (unloading) 

f  =  0,  A  =  0  (neutral)  (19) 

f  =  0,  A  >  0  (loading) 

to  cover  non-strain  hardening  materials. 

Finite  Element  Approach 

Various  types  of  axisyrnmetric  finite  elements  can  be 
derived  based  on  Equation  (1)  . 

2  3 

In  the  computer  program  described  in  the  Muraki 
Manual,  a  four  node-isoparametric  element  is  used.  Follow 
ing  a  standard  procedure  to  derive  an  element  stiffness 
matrix  and  element  forces,  the  coordinates  and  the  displac 
ments  for  the  elements  are  assumed  by 


(20) 


(21) 


38 

4 

v  =  X  h.v.  (21) 

i=l  1  1 

where  h^,  h^,  h^,  and  are  interpolation  functions,  x^ 

and  y.  are  nodal  coordinates,  and  u.  and  v.  are  nodal 
'i  ii 

displacements . 

A  key  of  the  above  derivation  is  use  of  the  same  in¬ 
terpolation  functions  for  both  the  element  coordinates  and 
the  element  displacements.  Details  of  the  interpolation 
functions  are 

hx  =  h  Cl  -  r)  (1  -  s) 

h_  »  h  (1  +  r)  (1  -  s) 

(22) 

h3  *  ^  Cl  +  r)  (1  +  S) 

h4  =  %(1  -  r)  (1  +  s) 

where  r  and  s  are  referred  to  as  natural  coordinates  and  va: 
from  -1  to  +1. 


The  element  strains  in  the  axisvmmetric  condition  are 
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Because  u  and  v  in  the  strain  are  a  function  of  x  and 
v,  the  following  relation  is  first  required  between  the  x-y 
and  r-s  systems 


f  \ 

3 

3x 

3r 

3r 

3r 

3 _ 

l  = 

3x 

5y 

3s 

N  / 

) 

3s 

3  s 

(24) 


Derivatives  of  the  displacements  with  regards  to  r  and 


s  are 


3u 

3r 


3u 

3s 


3u 

3r 


u 

s 


x  3h . 
y  _ i. 

1 1 

i=l3r 

LL  . 

i 

x  3h. 
j]  _ h. 

1 1 

i=i3r 

LL  . 

i 

x  3h . 

E  — — 
i=l3r 

v . 

i 

x  h . 

l 

i-1  r 

V  . 
1 

(25) 


Using  Equations  (24)  and  (25) ,  the  strains  in  Equation 
(23)  are  given  in  terms  of  nodal  displacements,  u^ ,  u2 ,  u3, 


and  u^  as  follows: 


{e}  =  | B  j { q } 


(26) 
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where 


r  •  1 T  r  •  1  , 

{t}  ~Ux'  “y '  lz'  -xy; 


xy 


{q}  =  {u1,  u2,  u3,  u4,  vv  v2,  v3,  v4) 


and  { B }  is  the  matrix  which  gives  a  relationship  between 
the  strains  and  the  ncdal  displacements.  It  is  noted  that 
{3}  is  a  function  of  r  and  s. 

The  stress-strain  relationship  of  Equation  (14)  is 
rewritten  in  a  simple  form  as: 

_  •  .  ,  I-*,.  _  •  0  _ 

{a}  =  D  { e }  +  {qw}  (27) 


where 


•xy 


} 


(D;  is  the  stress-strain  relationship  matrix  and  {r"}  is  the 
stress  vector  due  to  temperature  change,  and  the  temperature 
dependency  cf  material  properties  and  the  yield  function. 

The  element  stiffness  matrix  {k}  is  obtained  using  the 
principle  of  virtual  work  of  Equation  (1)  and  Castigliano '  s 
theorem  as, 

{k}  =  /.'/ {  3  } ^ { D}  ; 3 }v  d9dxdv  (28) 

2-  . '  \  3 } * { D } { B } y  j  J | drds 
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where 


3x 

3v 

3r 

3r 

3x 

ix 

3s 

3  s 

Introducing  a  numerical  integration,  we  have 


{k}  =  2-Z  {  3  }  ?  .  {  D }  {  3  }  .  .  I  J  I  .  .Y.  .a.  . 

i  a  13  i]  '  '  ID  13  ij 

^  f  J 


( 29 ) 


where  subscripts  i  and  j  indicate  that  Y^j ,  {3},  and  J !  are 

evaluated  at  coints  i  and  j.  cl.  .  is  weicht  factor  at  Doints 

13 

i  and  j  in  the  numerical  integration. 

The  program  in  the  Muraki  manual  uses  the  Gauss-Legendre 
numerical  integration  and  in  that  scheme  can  be  found  in 

text  books  on  numerical  integration. 

,  •  a 

The  element  force  vector  if  ;  corresponding  to  the  secon 
term  in  the  right  side  of  Equation  (27)  is  obtained  in  a 
similar  manner. 


!  f 


f  '  }  =  fff  {3} 1  {f  }yd9dxdv 

=  2r//{E}T{a"}yij| drds 


(30) 


Again  introducing  the  numerical  integration  we  have 


{T>  =  2-rZ  {3}T.{c5}y.  .  j  J;  .  .a.  . 

ifj  !3  13  13  13 


(31) 


In  the  current  version  of  the  computer  program  described 
in  the  Muraki  manual,  only  the  static  load  considered  is  the 
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is  the  concentrated  load  at  nodal  points.  Therefore  there 
is  no  other  element  force. 

The  stiffness  matrix  and  force  vector  for  the  whole 
body  (i.e.  the  whole  system)  are  obtained  by  assembling  the 
element  stiffness  matrices  and  force  vectors. 

The  equilibrium  equations  for  the  whole  system  are 
obtained  as 

Z{k] . {Q}  +  E{fe}.  =  Zif) .  (32) 

i  1  i  j  D 

where  i  and  j  stand  for  the  i-th  element  and  the  j-th  node 
respectively . 

{q}  is  the  vector  of  the  time  rate  of  nodal  displacemen 
{F;  is  the  vector  of  the  time  rate  of  concentrated  load. 

Equation  (32)  is  further  written  in  the  incremental 
form  as  follows: 


■  >  ’  : 
■ K ;  i  1 


:q; 


'  f  \p1 


(  *3 


3) 


where 

AQ  =  QAt 

A  £  ^  1  \  4_ 

AF  =  FAt 

At  is  the  time  increment. 

It  is  noted  that  the  time  increment  is  included  in 
Equation  (33)  implicitly. 
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Once  the  nodal  displacements  for  the  system,  are 

obtained,  the  element  strains  are  calculated  by  Equation  (26) 
and  then  the  corresponding  stresses  by  Equation  (27) . 

>7hen  there  is  no  further  increment  of  concentrated  and/ 
or  thermal  loads,  the  analysis  is  complete. 

When  there  is  another  increment  of  the  loads ,  new 
equilibrium  equations  are  calculated  using  the  nodal  dis¬ 
placements  and  stresses  obtained  in  the  preceding  calculation. 
This  procedure  of  incremental  analysis  continues  until  the 
final  increment  is  solved. 


m 


m 
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2 .  Geometry  of  Weld  Specimen  and  Finite  Element 
Representation  of  Weldment  used  by  Muraki  Program: 

Finite  element  representation  of  the  weldment  is 
depicted  in  Figure  V-2.  A  comparison  of  Figure  V-2  with 
the  finite  element  representation  of  the  weldment,  Figure  V-l, 
used  in  the  heat  flow  program  will  shew  that  an  axis  trans¬ 
formation  as  well  as  a  shift  of  the  origin  were  required  in 
order  to  use  heat  flow  program  output  as  input  for  the  Muraki 
Program. 

For  the  Muraki  program,  the  v-axis  corresponds  to  the 
radial  direction  and  the  x-axis  is  along  the  cylinder  length. 
The  z-axis  corresponds  to  the  direction  of  welding  arc  travel, 


normal  to 

the 

x-y 

plane  as 

shown  in  Figure  V-2,  pcsi 

tive 

direction 

out 

of 

the  paper . 

The  origin  (x=y=0)  for 

the  Muraki 

program  analysis 

correspond 

s  to  a  point  at  the  cents 

r  cf  the 

cylinder  end.  2=0  corresponds  to  the  time  cf  arc  initiation. 

The  Muraki  program  also  required  the  specif i ration  cf 
nodal  coordinates  (x,y  values) .  Element  connectivity  was  input 
to  the  Muraki  program  by  identification  cf  the  nodes,  in  counter¬ 
clockwise  fashion  beginning  with  the  node  closest  to  the  heat 
source,  associated  with  each  element.  Medal  coordinates  input 
to  the  Muraki  program  are  listed  in  Table  V-3  and  element  con¬ 
nectivity  is  shown  in  Table  V-C . 

Therncphysical  material  properties  required  by  the 
Muraki  program  were  taken  from  Figures  II- 1  throuen  XI-lv 
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in  Chapter  II  of  this  study.  Welding  parameters  required  as 
input  data  for  the  first  pass  of  the  girth  weld  were  the  same 
as  those  used  by  the  heat  flow  program.  These  parameters 
were  discussed  and  values  presented  earlier  in  this  chapter. 


Y  (iKUlfcb)  FINITE  ELEMENT  REPRESENTATION  OF  WELDMENT  GEOMETRY 

FOR  STRAIN,  STRESS,  AND  DISTORTION  ANALYSIS 


FIGURE 


NODE  COORDINATES  (MURAKI  PROGRAM) 


n 
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CHAPTER  VI 


Comparison  of  temperatures  and  values  of  strain,  stress, 
and  distortion  during  production  of  the  HY-130  cylinder  girth 
weld  with  those  predicted  by  use  of  the  heat  flow  and  Muraki 
programs  follows  in  the  form  of  tabulated  data  and  in  graph¬ 
ical  presentations.  During  production  of  the  girth  weld,  arc 
travel  began  at  the  seam  weld  opposite  the  cylinder  seam  about 
which  thermocouples  and  strain  gages  were  placed.  Times  at 
which  data  comparisons  are  presented  were  selected  as  repre¬ 
senting  a  thermal  cycle  during  which  temperature,  strains, 
stress,  and  distortion  underwent  greatest  changes. 

Temperatures  determined  experimentally  at  the  eight 
thermocouple  locations  for  the  times  of  interest  are  shown 
in  Table  VI-A.  Distances  from  weld  centerline  are  indicated 
beside  the  column  headings  designated  Tl,  T2,  T3,  T4,  T4', 

T31,  T2 ' ,  and  Tl'.  These  temperatures  are  presented  graph¬ 
ically  in  Figure  VI-1.  Temperature  distribution  predicted  by 
use  of  the  heat  flow  program  is  presented  by  a  computer  gener¬ 
ated  graph,  Figure  VI-2. 

Qualitative  comparison  of  the  analytical  and  measured 
temperature  distributions  shows  remarkable  coincidence.  The 
slightly  higher  overall  values  of  the  analytically  generated 
values  may  be  explained  by  the  assumed  welding  arc  efficiency 
(narc)  of 
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FIGURE  VI-1  GRAPH  OF  EXPERIMENTALLY  DETERMINED  TEMPERATURES 
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For  the  purpose  of  presenting  comparative  strain,  stress, 
and  distortion  data,  the  following  terms  are  defined: 

STR ( * . * ) X  =  EXPERIMENTALLY  DETERMINED  STRAIN 

(MICROSTRAIN)  AT  A  POINT  *.*  INCHES 
FROM  THE  WELD  CENTERLINE. 

STR ( * . * ) A  =  ANALYTICALLY  PREDICTED  STRAIN 

(MICROSTRAIN)  AT  A  POINT  *.*  INCHES 
FROM  THE  WELD  CENTERLINE. 

STS  (.*.*)  X  =  EXPERIMENTALLY  DETERMINED  STRESS  (lbf/.  2) 

in 

AT  A  POINT  *.*  INCHES  FROM  THE  WELD 
CENTERLINE . 

STS  (*  .  *)  A  =  ANALYTICALLY  PREDICTED  STRESS  (lbf/^2) 

AT  A  POINT  *.*  INCHES  FROM  THE  WELD 
CENTERLINE . 

DST ( * . * ) X  =  EXPERIMENTALLY  MEASURED  RADIAL  DISTORTION 
IN  INCHES  AT  A  POINT  *.*  INCHES  FROM  THE 
WELD  CENTERLINE. 

DST ( * . * ) A  =  ANALYTICALLY  PREDICTED  RADIAL  DISTORTION 
IN  INCHES  AT  A  POINT  *.*  INCHES  FROM  THE 
WELD  CENTERLINE. 

Positive  values  of  strain  and  stress  correspond  to  the 
tensile  condition  and  negative  values  correspond  to  compressive 
strains  and  stresses.  Positive  values  of  radial  distortion 
correspond  to  displacement  in  an  outward  direction  from  the 
cylinder  center.  Tables  and  graphs  on  the  following  pages 
use  the  terns  as  defined  abo^e  to  compare  data.  Strain  and 
Stress  directions  correspond  to  directions  shewn  in  Figure  V-2. 
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Tables  VI-B  through  VI-E  present  in  tabular  form 
comparisons  between  strains  measured  during  production  of 
the  HY-130  cylinder  girth  weld  and  those  predicted  by  use 
of  the  Muraki  program.  Tables  VI-F  through  VI-I  present 
the  analytically  predicted  stresses.  Experimentally  deter¬ 
mined  stresses  were  not  tabulated  since  experimentally 
determined  strains  could  not  be  transformed  into  stresses. 

There  are  three  (ignoring  creep)  components  of  strain;  elastic, 
plastic,  and  thermal.  Elastic  and  thermal  constitutive  rela¬ 
tions  are  simple  and  the  stresses  resulting  therefrom  could 
be  calculated  if  these  strains  could  be  separated  into  their 
components.  However,  the  constitutive  relationships  among 
elastic,  plastic,  and  thermal  strains  are  complex,  made  further 
so  by  the  presence  of  a  strain  hardening  material  and  by  the 
dependence  of  all  material  properties  on  temperature.  For 
this  reason,  experimentally  determined  strains  could  not  be 
transformed  into  stresses.  Values  of  radial  distortion  measured 
during  production  of  the  girth  weld  and  predicted  by  the  Muraki 
program  are  shown  in  Table  VI -J.  Comparative  graphs  of  measured 
and  predicted  data  are  shown  in  Figures  VI-3  through  VI-9. 


LONGITUDINAL  STRAIN 
(X-DIRECTION) 


TRANSVERSE  STRAIN 
(Z -DIRECTION) 


TIME (sec) 


125.0 


127.0 


STR ( 1 . 0 ) X  |  STR ( 1 . 0 ) A 


STR ( 1 . 0 ) X  I  STR (1.0) A 


145.0 


150.0 


155.0 


160.0 


170.0 


200.0 


TABLE  VI-3  EXPERIMENTAL  AND  ANALYTICAL  STRAIN  1.0  IN  FROM. 
WELD  CENTERLINE 
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LONGITUDINAL  STRAIN 
(X-DIRECTION) 

TRANSVERSE  STRAIN 
(Z-DIRECTION 

TIME (sec) 

STR ( 1 . 5 ) X 

STR (1. 5) A 

STR (1. 5) X 

STR (1.5) A 

125.0 

-245 

0 

320 

0 

127.0 

-310 

0 

365 

0 

123.0 

-325 

0 

395 

1 

129.0 

-345 

2 

410 

16 

129.5 

-345 

-70 

425 

27 

130.0 

-345 

-21 

450 

41 

130.5 

-340 

-105 

470 

92 

131.0 

-340 

-259 

485 

185 

132.0 

-330 

-358 

505 

260 

133.0 

-330 

-318 

570 

263 

134.0 

-320 

-297 

626 

277 

135.0 

-325 

-288 

681 

293 

136.0 

-330 

-283 

706 

311 

137.5 

-335 

-265 

750 

336 

140.0 

-340 

-205 

749 

373 

145.0 

-330 

11 

664 

433 

i 

150.0 

-350 

276 

560 

477 

155.0 

-400 

527 

496 

513 

160.0 

-390 

651 

427 

530 

170.0 

-330 

1060 

257 

590  | 

200.0 

-291 

1490 

-246 

- ! 

666  j 

TABLE  VI-C  EXPERIMENTAL  AND  ANALYTICAL  STRAIN  1.5  IN  FROM 
WELD  CENTERLINE 


9 


10 


LONGITUDINAL  STRAIN 
(X-DIRECTION) 

TRANSVERSE  STRAIN 
( Z -DIRECTION) 

TIME (sec) 

STR ( 2 . 0 ) X 

STR (2 . 0) A 

STR ( 2 . 0 ) X 

STR (2.0) A 

125.0 

-310 

0 

315 

0 

127.0 

-355 

0 

385 

0 

128.0 

-360 

0 

395 

1 

129.0 

-370 

-5 

410 

14 

129.5 

-385 

-16 

435  . 

22 

130.0 

-400 

-32 

440 

33 

130.5 

-395 

-115 

445 

68 

131.0 

-390 

-264 

455 

132 

132.0 

-380 

-365 

460 

186 

133.0 

-365 

-334 

475 

193 

134.0 

-355 

-321 

490 

207  1 

135.0 

-320 

-320 

498 

221  j 

136.0 

-310 

-325 

500 

236 

137.5 

-283 

-331 

502 

257  | 

1 

140.0 

-225 

-337 

520 

288 

145.0 

-220 

-304 

555 

338  | 

150.0 

-384 

-216 

606 

378 

155.0 

-509 

-94 

666 

413 

160.0 

-636 

-16 

689 

431 

170.0 

-625 

304 

600 

496  | 

200.0 

-326 

836 

164 

587 

TABLE  VI-D  EXPERIMENTAL  AND  ANALYTICAL  STRAIN  2.0  IN 
FROM  WELD  CENTERLINE 
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LONGITUDINAL  STRAIN 
(X-DIRECTION) 

TRANSVERSE  STRAIN 
( Z-DIRECTION) 

TIME (sec) 

STR (2.5) X 

STR (2. 5) A 

STR ( 2 . 5)  X 

STR ( 2 . 5)  A 

125.0 

-175 

0 

295 

0 

127.0 

-220 

0 

310 

0 

128.0 

-250 

0 

320 

1 

129.0 

-265 

-5 

325 

*11 

129.5 

-265 

-8 

335 

16 

130.0 

-265 

-13 

355 

24 

130.5 

-265 

-30 

355 

47 

131.0 

-265 

-61 

365 

88 

132.0 

-265 

-86 

365 

124 

133.0 

-260- 

-85 

365 

132 

134.0 

-260 

-89 

355 

143 

135.0 

-255 

-93 

355 

154 

136.0 

-245 

-98 

350 

166 

137.5 

-223 

-104 

340 

'  '  i 

182 

140.0 

-125 

-114 

340 

205 

) 

145.0 

-us 

-123 

355 

244 

150.0 

-249 

-118 

457 

276 

155.0 

-379 

-99 

552 

305 

160.0 

-501 

-83 

625 

320  | 

170.0 

-560 

1 

624 

377 

200.0 

-376 

203 

365 

463  | 

TABLE  VI-E 


EXPERIMENTAL  AND  ANALYTICAL  STRAIN  2, 
FROM  WELD  CENTERLINE 


5  IN 


j 


LONGITUDINAL  STRESS 

TRANSVERSE  STRESS 

(X-DIRECTION) 

(Z-DIRECTION) 

TIME ( sec)  !  STS ( 1 . 0 )  X  |  STS (1.0) A  |  STS (1.0) X  |  STS (1.0) A 


125.0 


127.0 


123.0 


129.0 


129.5 


130.0 


130.5 


131.0 


137.0 


133.0 


134.0 


135.0 


136.0 


137.5 


140.0 


5 . 0 


150.0 


155.0 


160.0 


170.0 


200.0 


0 


3 


143 


436 


207 


-295 


■2658 


■6614 


-9464 


-8670 


■8263 


-8235 


■8420 


-8619 


-8727 


-7337 


-4913 


-2501 


-1321 


2508 


6329 


0 


22 


374 


593 


888 


1235 


2669 


5369 


11431 


6405 


5400 


3860 


1908 


-1588 


-7665 


■16930 


■21859 


■24190 


■24732 


•25156 


■21790 


TABLE  VI-F  PREDICTED  X  AND  Z  DIRECTION  STRESSES  1.0  IN 
FROM  WELD  CENTERLINE 


TABLE  VI-G  PREDICTED  X  AND  Z  DIRECTION  STRESSES  1.5  IN 
FROM  WELD  CENTERLINE 
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LONGITUDINAL  STRESS 
(X-DIRECTION) 


TRANSVERSE  STRESS 
(Z -DIRECTION) 


-3  87  2 


TABLE  VI-H  PREDICTED  X  AMD  Z  DIRECTION  STRESSES  2.0  IN 


LONGITUDINAL  STRESS 

TRANSVERSE  STRESS 

(X-DIRECTION) 

( Z -DIRECT ION) 

STS (2. 5) X  STS (2. 5} A 


0 


0 


-2 


-23 


-50 


-89 


-269 


-595 


-827 


-780 


-772 


-782 


-807 


-842 


-889 


-958 


-1025 


-1140 


-1231 


-1689 


-2710 


STS (2. 5) X  |  STS  (2. 5) A 


17 


295 


450 


656 


1301 


2465 


3487 


3687 


4001 


4299 


4618 


5036 


5658 


6567 


7051 


7139 


7013 


5963 


2510 


TABLE  VI-I  PREDICTED  X  AMD  2  DIRECTION  STRESSES  2,5  IN 
FROM  WELD  CENTERLINE 


RADIAL  DISTORTION 
(Y-DIRECTION) 


RADIAL  DISTORTION 
(Y-DIRECTION) 


TIME (sec) 


125.0 


127.0 


128.0 


129.0 


129.5 


130.0 


130.5 


131.0 


133.0 


134.0 


135.0 


136.0 


137.5 


0. 


145.0 


150.0 


155.0 


160.0 


170.0 


200.0 


DST ( 0 . 5 ) X  I  DST ( 0 . 5 ) A  DST(1.0)X  IDST(1.0)A 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.001 


0.001 


0.001 


0.002 


0.002 


0.003 


0.003 


0.003 


0.003 


0.003 


0.002 


0.002 


0.0000 


0.0000 


0.0000 


0.0000 


0.0002 


0.0002 


0.0009 


0.0026 


0.0035 


0.0029 


0.0028 


0.0028 


0.0029 


0.0032 


0.0035 


0.0039 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0041 

0.004 

0.0043 

0.004 

0.0044 

0.004 

0.0046 

0.004 

0.0049 

0.003 

•DEFLECTION  TOO  SMALL  TO  QUANTIFY 

TABLE  VI-J  EXPERIMENTAL  AND  ANALYTICAL 
1.0  IN  FROM  WELD  CENTERLINE 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0005 


0.0013 


0.0028 


0.0038 


0.0037 


0.0037 


0.0037 


0.0038 


0.0039 


0.0041 


0.0044 


0.0045 


0.0046 


0.0047 


0.0049 


0.0052 


DISTORTION  0.5  IN  AND 


0.001 


0.001 


0.001 


0.001 


0.002 


0.002 


0.003 


0.004 


FIGURE  VI-3  COMPARATIVE  PLOT  OF  EXPERIMENTAL  AND  ANALYTICAL  STRAIN  1.0  IN  FROM 
WELD  CENTERLINE 
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FIGURE  VI-9  COMPARATIVE  PLOT  OF  EXPERIMENTAL  AND  ANALYTICAL  RADIAL  DISTORTION 
0.5  IN  AND  1,0  IN  FROM  WELD  CENTERLINE 
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CHAPTER  VII 

CONCLUSIONS 

This  study  of  residual  stresses  and  distortion  in  high, 
strength  steel  weldments  has  shown  that  existing  computer 
programs  require  further  refinement  in  order  to  be  used  as  a 
day-to-day  tool  of  the  welding  engineer. 

The  occurrence  of  residual  stresses  and  distortion  during 
production  of  a  weldment  is  a  complex  phenomena.  Past  efforts 
to  simulate  the  welding  process  involved  clever  simplifying 
assumptions  followed  by  attempts  to  carry  out  experiments  which 
supported  these  assumptions.  The  new  approach  to  welding 
simulation  employs  the  enormous  power  of  computer  techniques 
to  analyze  what  actually  takes  place  during  the  welding  process. 
However,  if  a  large  computational  process  is  used  to  model  the 
complex  welding  process,  the  thermophysical  material  properties 
on  which  the  analysis  is  based  must  be  accurately  known. 

The  axisymmetric  finite  element  program  which  was  developed 
by  Muraki  for  analysis  of  thermal  stress  and  metal  movement 
during  welding  is  promising  in  that  it  showed  good  qualitative 
agreement  with  carefully  generated  experimental  data.  Heat 
flow  and  resulting  temperature  distributions  within  weldments 
can  be  accurately  predicted  by  existing  methods.  Qualitative 
and  quantitative  comparison  between  analytical  and  measured 
temperature  distributions  showed  remarkable  coincidence. 
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CHAPTER  VIII 
RECOMMENDATIONS 

Inasmuch  as  analytical  methods  to  predict  stresses  and 
distortion  resulting  from  thermal  cycles  can  only  be  as  good 
as  the  thermophysical  material  properties  on  which  they  are 
based,  it  is  recommended  that  studies  be  undertaken  to 
accurately  determine  the  temperature  dependence  of  high- 
strength  steel  properties  at  alevated  temperatures.  This 
would  enhance  further  efforts  to  simulate  the  welding 
process , 

For  further  study  of  computer  simulation  of  multi-pass 
welding,  it  is  recommended  that  a  micro-processor  be  used 
to  process  welding  parameter  data.  Reduction  of  strains 
requiring  temperature  compensation  is  a  tedious  task, 
especially  if  measurements  are  made  at  several  locations 
and  recorded  on  one  photo  tape.  This  additional  source  of 
error  could  be  reduced  or  eliminated  if  parameters  of 
interest  could  be  calculated  internally  and  displayed 
directly. 

High  temperature  strain  gages  could  significantly  con¬ 
tribute  toward  a  more  meaningful  study.  The  400  °F  limitation 
of  gages  used  in  this  study  severely  restricted  '  their  loca¬ 
tion. 
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